Percolation in Models of Thin Film Depositions 
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We have studied the percolation behaviour of deposits for different (2+l)-dimensional models of 
surface layer formation. The mixed model of deposition was used, where particles were deposited 
selectively according to the random (RD) and ballistic (BD) deposition rules. In the mixed one- 
component models with deposition of only conducting particles, the mean height of the percolation 
layer (measured in monolayers) grows continuously from 0.89832 for the pure RD model to 2.605 
for the pure RD model, but the percolation transition belong to the same universality class, as 
in the 2- dimensional random percolation problem. In two- component models with deposition of 
conducting and isolating particles, the percolation layer height approaches infinity as concentration 
of the isolating particles becomes higher than some critical value. The crossover from 2d to 3d 
percolation was observed with increase of the percolation layer height. 

PACS numbers: 64.60.Ak, 72.60.-fg, 81.15.Aa, 89. 75. Da 



The thin film formation processes by deposition of par- 
ticles on a substrate are of great interest both from the- 
oretical as well as experimental point of view ^ . The 
different aspects of this problem are important in the 
technical applications for production of thin-film devices, 
metal-insulator mixture films, composite films with spe- 
cific physical properties, etc. 

The rather important field of investigation is related 
to the electrical conductivity of thin films, which de- 
pends strongly on their morphology and microstructure. 
Many works were devoted to investigations of the frac- 
tal, percolation and electrical properties of thin films and 
deposits It was shown that the percolation transi- 
tion in very thin (quasi 2-dimensional) films belongs to 
the same universality class as in the random percolation 
problem l^j, ^. The film electrical conductivity shows 
also a clear transition from the two-dimensional to three- 
dimensional behaviour with film thickness increase ^ . 
Some correlations were observed between the conductiv- 
ity and porosity for deposits grown in a model of ballistic 
deposition . Jensen et al. and Family [jll] investi- 
gated in their numerical simulation works the percolation 
behavior for different models of submononolayer deposits 
on two-dimensional substrates. 

The purpose of this work is to study the percolation be- 
haviour for different lattice models of three dimensional 
deposits growing on the plane substrates. The spanning 
cluster forms in the substrate plane. The percolation in 
deposits has a correlated character, because the sites of 
lattice get filled dynamically during the growth in accor- 
dance with the deposition rules. 

In our model the particles are modeled by unit cubes. 
They are deposited on an initially flat horizontal surface 
on the x — y plane of size Lx L. The particles come down 
vertically along the ~z direction with the integer x, y co- 
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FIG. 1: Scheme of percolation cluster formation for 2-1-1 
dimensional deposition model. 



ordinates and are deposited on the substrate either by 
the ballistic deposition (BD) process or by the random 
deposition (RD) process or by a mixture of both the pro- 
cesses. In the RD process a particle comes down verti- 
cally till it lands over a particle on the substrate where 
as in the BD process the particle gets stuck to the sub- 
strate when any of its four vertical sides comes in contact 
with any previously deposited particle of the substrate or 
it directly lands on the substrate. In a mixed RD-I-BD 
model, the s fraction of particles is deposited according 
to the BD-rule and the remaining 1 — s fraction is de- 
posited according to the RD-rules. We call this model 
as: BDsRDi_s model. The parameter s allows a contin- 
uous tuning between RD model with no interaction be- 
tween particles and BD model with strong short-ranged 
interaction between particles. The model of this type or 
other similar models where different kind of interactions 
between particles may exist are widely used for simula- 
tion of structure of thin films with realistic morphology 
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FIG. 2: Plots of height hoo and density poo of deposit in per- 
colation point versus parameter s for one-component mixed 
BDsRDi-s model. The data error is of order of data symbol 
size. The lines serve as a guide to the eye. 



In the mono-component model all particles are consid- 
ered to be conducting. In the bi-component model we 
have a / fraction of insulating particles and the 1 — / 
fraction of conducting particles. 

Particles are deposited on the substrate one after an- 
other and the average height of the deposit grows. Con- 
duction takes place between two conducting particles 
when they have one surface in contact. We stop the 
growth process when the deposit starts conducting for 
the first time in the direction parallel to the surface. At 
this percolation point, a spanning cluster across the sys- 
tem is formed along the x or y direction (Fig. 0).The per- 
colation point is easily checked by a Hoshen-Kopelman 
algorithm ||l3|] . 

During the deposition process, the time elapsed is mea- 
sured in units of the number of equivalent complete layers 
deposited. Therefore N particles have been deposited in 
time t = N/L'^. On the other hand the mean height of the 
deposit at time t is. h — Yix,yhxy{i) / L"^ . The percolation 
density is the volume fraction p of the conducting parti- 
cles at the percolation point i.e., the ratio of the number 
of conducting particles and the total volume of the de- 
posit p = Nc/ {JiLP'). In mono-component model Nc = N 
and p = t/h. For RD model the bulk of the deposit 
is compact (without any pores in the vertical columns) 
but it has a rough interface. Therefore the RD limit at 
s = corresponds to p — 1. For ballistic deposition (BD) 
model at s = 1 the deposit is porous and consequently 
p < 1. For a bi-component model we have a mixture of 
conducting and insulating particles. Here the total den- 
sity of the deposit including conducting and insulating 
particles is ptotai — N/ihL'^) and Nc — N{1 — /) and 

p = t{l- f)/h^Ptotal{l- f). 

For finite extensions [L) of the substrates the percola- 
tion height h{L) and p{L) are L dependent. The values 
of h{L) and p{L) are determined for different substrate 



sizes L varied from 8 to 2048 and the periodical bound- 
ary conditions were applied in deposition rules along the 
directions x and y. Results were averaged over 100-5000 
different runs, depending on the size of the lattice and 
required precision. 

In analogy with the corresponding finite size be- 
haviours in the ordinary percolation, we assume the fol- 
lowing relations: 



p{L) = pac. + apL 



-i/i/p 



and 



(1) 



(2) 



where Up, Vh ~ 4/3. 

The probability that a particle is deposited along a 
particular vertical line on the x — y plane is 1 /L^ which 
is very small when L is large such that the mean height 
h is maintained at a fixed value when TV particles are 
deposited. This implies that in the RD model the number 
h of particles in an arbitrary column of particles follow a 
Poisson distribution: 



P{h) = {e-h{h'')/hl. 



(3) 



Therefore the probability of an empty column (h = 0) is 
equal to P(0) = e^'*. In percolation point 



PiQ)^l-p2d^e-'\ 



(4) 



where p2d is the percolation threshold for the square lat- 
tice site percolation problem. 

Taking into account the finite size scaling behavior of 

p: 



P2d{L) = P2d,oo + a2d,pL ^/"p. 



(5) 



where p2d.oo = 0.592746.. is the percolation concentration 
in the limit of infinite system (L — > oo) and i/p = 4/3 is 
a correlation length scaling exponent |14| , we obtain 



h{L) = - ln(l - p2d,oo - a2d,pL ^Z"") 



(6) 



= /loo +ln(l - 



Q2d,p 



{1 - P2d,oo)L^/'''' 



■ah 



L-^/''\ (7) 



where hoo = — ln(l 

P2d,oo), Vh = Vp- 



P2d.. 



0.89832, ah = a2d,p/(l 



We see that for pure RD model Vh 



4/3. 



So, the RD model belong to same class of universality 
as the 2-dimensional random percolation model. This 
fact reflect the small mean height of (2+l)-dimensional 
random deposit hoc ~ 0.89832, which is only slightly 
higher than mean height of 2-dimensional random de- 
posit hoc « 0.592746. 

On the basis of numerical simulations we estimate that 



for pure BD at s = 1, /lo 



2.605 ± 0.005 and Po 
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FIG. 3: Plots of height hoc and density of deposit in 
percolation point versus fraction of isolating particles for two- 
component BD and RD models. In cases when it is not show 
directly the data error is of order of data symbol size. The 
solid lines serve as a guide to the eye. Vertical dashed lines 
show the critical concentrations which are /c = 0.227 ± 0.001 
and /c = 0.70 ± 0.01 for BD and RD models respectively. 
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FIG. 4: Plots of scaling exponents Vh and in Eqs. 
versus fraction of isolating particles for two-component BD 
and RD models. The solid lines serve as a guide to the eye. 
Vertical dashed lines show the critical concentration which 
are /c = 0.227 ± 0.001 for BD model and fc = 0.70 ± 0.01 for 
BD model. 



0.620±0.005 in the limit L = oo. Using these asymptotic 
values we plot hoo — h{L) and poo —p{L) vs, L on double 
logarithmic scales and in both cases plots correspond to 
the slope of -3/4, which means Vh = Vp — 4/3. 

Percolation height hoo and the percolation densities Poo 
are similarly calculated for the mixed BDsRDi„s model 
varying the mixing parameter s and plotted in Fig. |^. 
The height of the deposit hoo increases and its density 
decreases smoothly with increase of the fraction of de- 
posited BD particles. In the limit of pure RD model, the 
theoretical value hoo ~ 0.89832 is observed well. Same as 
both pure BD and RD models, the mixed BDsRDi_s also 
displays the scaling behaviour described by Eqs. ([l]j^) 
with scaling exponent Vh = Vp ~ 4/3. Using the cal- 



FIG. 5: Log-log plot of pao — Pc,oo versus hoo for two- 
component BD(squares) and RD(triangles) models. See the 
text for the details. 



culated hoo{s) and Poo(s) dependencies and substituting 
Vh = Vp = 4/3 into Eqs. (flS the coefficients au and ap 



versus s were obtained. Both of these coefficients au and 
ap increase with s. It is important to note that for the 
pure RD model = 0. It means that there is no finite 
size scaling for the RD model, as for any L at p = 1 the 
deposit is compact without pores by definition. 

We can conclude that the mixed BD^RDi-s model, 
presumably, belong to same class of universality as the 
2-dimensional random percolation model at any value of 
s. 

Figure || presents the deposit height hoo and density 
of conducting particles Poo in the percolation point esti- 
mated in the limit of L — > cx) versus fraction of isolating 
particles /. The dependencies of hoo{f) show the typical 
percolation behaviour: as / reaches some critical value 
/c the value of hoo goes to infinity; it means that there 
is no percolation at any finite height of deposit. The es- 
timated values of critical concentrations of the isolating 
particles arc fc{BD) = 0.227 ± 0.001 for the BD model 
and fc{RD) = 0.70 ± 0.01 for RD model. 

For the RD model, the total density of particles is 
Ptotai = 1 by definition, and, so, the density of con- 
ducting particles is p = ptotai{^ - f) = 1 - /• The 
linear law of Poo decrease with / increase is actually ob- 
served in simulation data for the RD model, this law is 
rather close to linear for the BD model (Fig. ||). In the 
critical point / = /c, the density of conducting parti- 
cles is equal to Pc,oo = 0.232 ± 0.001 for the BD model, 
and Pc,cx3 = 0.30 ± 0.01 for the RD model. This value 
for the RD model is very close to percolation concentra- 
tion for the random percolation on a simple cubic lattice 
p = 0.311609 ||l^. The estimated value of the total den- 
sity of deposit ptotai.,oo = 0.300 ±0.001 for the BD model 
coincides with the previously reported data for the de- 
posit density extrapolated to the infinite-system limit for 
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the BD model |6). 

The scaUng exponents Vh, obtained by the least 
square fit of Eqs. (|l|,||) versus / are presented at Fig. 
^ for the BD and RD models. All fitting procedures were 
done at the fixed values oipoo, two free parameters and 
correlation coefficients were higher than 0.998. For the 
RD model, the value of p is independent from system size 
L by definition, so, only dependency is presented at 
Fig.|. 

The scaling exponent Vh for the BD model contin- 
uously grows with s from Vh ~ 4/3 characteristic for 
single- component model to Vh ~ 2.4 near critical point 
f ^ fc = 0.227. The dependence of Vh versus / 
for the RD model is less pronounced. The correlation 
length exponent j/p for BD model smoothly decreases 
from Vp = 1/2(1 = 4/3 characteristic for the 2d-systems 
(at f = 0) to Vp = v^d — 9/10 characteristic for the 
3d-systems (at f = fc = 0.227) This behaviour 

can be easily understood, as far as the height of per- 
colation deposit at high / increases and approaches the 
infinity at the critical concentration of isolating particles. 
So, the BD-model percolation of deposit in the vicinity 
of / = /c = 0.227 may be considered as a percolation 
in a three-dimensional system and the universality class 
of this percolation model is presumably the same as for 
the random percolation model. Existence of the 2d-to- 
3d percolation crossover was experimentally observed in 
random metal-insulator mixture films, when thickness of 
the films deposited was increased [Q. 

As far as at / < /c, {L ^ oo) the height of a perco- 
lation cluster hoc remains finite even for infinitely large 
substrate dimension, it is interesting to check poo versus 
/loo for existence of scaling: 

Poo = Pc,oo + apji^^""'" ■ (8) 

Fi gure ^ shows the log- log presentation of p^o ~ Pcoo 
versus versus hoo for the BD and RD models. For the BD 
model, we put Pc,oo = 0.232. The solid line 1 corresponds 
to best fit of Eq. (||) to the data (filled squares) with pa- 
rameters Vph = 1.20 ± 0.01 and aph = 0.883 ± 0.001. 

Putting the value Pc.oo — 0.311609 for the RD model, 
which is equal exactly to the concentration for the ran- 
dom percolation problem on simple cubic lattice, we ob- 
tain data represented by the open triangles. The dashed 
line 2 is drawn as a guide to the eyes and the scaling is 
rather poor. But when we put a somewhat higher value 
Pc,oo = 0.30 (remind that the value obtained from h^o 
versus / dependencies is Pc,oo = 0.30 ± 0.01) we get the 
data represented by filled triangles. Now the scaling is 
rather good and the solid line 3 corresponds to best fit 
of Eq. (H) to the data (filled triangles) with parameters 
Vph = 1.04 ± 0.01 and Qpu = 0.633 ± 0.004. 

In summary, we have investigated the percolation in 
the direction parallel to the surface for different models 
of deposition layer formation. In one-component models 



with deposition of only conducting particles, the height 
of deposits in the point of percolation is finite and is 
hoo = 0.592746 for the RD model hoo = 2.605 for the 
BD model. The mixed BDsRDi_s model, presumably, 
belongs to the universality class of 2d-random percolation 
problem. In two-component models with conducting and 
isolating particles, the percolation layer height hoo can 
be varied in wide range by tuning of a concentration of 
the isolating particles / and a clear crossover from 2d 
to 3d percolation is observed with increase of hoo- The 
percolation in layer is impossible when / exceeds some 
critical concentration fc- The weakening of inter-particle 
interaction results in increasing of threshold value of /c, 
from 0.227 for BD model « 0.7 for RD model. 
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